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ABSTRACT 

The tidal breakup of binary star systems by the supermassive black hole (SMBH) in the center of the 
galaxy has been suggested as the source of both the observed sample of hypervelocity stars (HVSs) in 
the halo of the Galaxy and the S-stars that remain in tight orbits around Sgr A*. Here, we use a post- 
Newtonian N-body code to study the dynamics of main-sequence binaries on highly elliptical bound 
orbits whose periapses lie close to the SMBH, determining the properties of ejected and bound stars 
as well as collision products. Unlike previous studies, we follow binaries that remain bound for several 
revolutions around the SMBH, finding that in the case of relatively large periapses and highly inclined 
binaries the Kozai resonance can lead to large periodic oscillations in the internal binary eccentricity 
and inclination. Collisions and mergers of the binary elements are found to increase significantly for 
multiple orbits around the SMBH, while HVSs are primarily produced during a binary's first passage. 
This process can lead to stellar coalescence and eventually serve as an important source of young stars 
at the Galactic center. 

Subject headings: black hole physics-Galaxy:center-Galaxy:kinematics and dynamics-stellar dynamics 



1. INTRODUCTION 

The center of the Milky Way can be regarded as a 
rather unusual galactic nucleus in many respects. The 
supermassive black hole (SMBH) at the Galactic center 
has a mass of Mmw 4 X IO^Mq (|Ghez et al. 1 120081 : 
iGillessen et al.l |2009() , making it per haps the smallest 
SMB H with a well-determined mass (|Ferrarese fc FordI 
I2005f ). Furthermore, while most known galactic nuclei 
have relaxation times much longer than the Hubble time, 
the relatively high density at the center of the Milky Way 
( ^ 10^ stars pc~'^) and the small mass of the SMBH 
imply a two-body relaxation time that is a s short as a 
few G yr inside the SMBH influence radius (Alexander I 
120051 ). The combination of high stellar densities, large 
velocities, and short relaxation time effectively make the 
inner parsec of the Milky Way a coUisional system. 

In a coUisional environment, otherwise rare dynam- 
ical processes can take place at an appreciable rate. 
These include dynamical encounters between stars, bi- 
naries and higher-order systems. If a massive black 
hole is present, as is believed to be the case for most, 
if not all, galactic nuclei, stars and binaries can in- 
teract with it and be ejected with extreme veloci- 
ties. Gravitational encounters involving stellar bina- 
ries and the SMBH have recently been studied by 
several a uthors in the context o f hypervelocity stars 
(HVSs) (|Gualandris et all [20051 : IBromlev et all [2005 
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IGinsburg fc Loebl l200l ISesana et al.1 l2007i) . HVSs 
have such extreme velocities (500 — 1000 kms~^ as mea- 
sured in the Galactic halo) that they require a dynam- 
ical encou nter with a massive black h ole for their ex- 
planation ([Hills lllQSSHYu fc Tremgdliillfoda ). The idea 
is that stellar binaries on low angular momentum or- 
bits interact with the SMBH and are dissociated by its 
strong tidal field. As a result, one component is cap- 
tured by the SMBH into a wide eccentric orbit while the 
other star is ejected to infinity with large velocity. After 
the f ortunate discovery of the first HVS in the Galactic 
halo (jl rown et ani2005n . a handful have been found es - 
caping the Galaxy at large speeds ([Brown et al.ll2006l ). 
More recent surveys of hypervelocity stars found 16 stars 
with veloci ties larger than the Galactic escape velocity 
([Brown et al. 2007, 2009). 

Given the high stellar density at the Galactic center, 
it is plausible that stars interacting with the SMBH will 
come close enough to each other that finite size effects 
become important and physical collisions occur. In this 
case, the SMBH acts a s catalyst for their interaction. 
IGinsburg fc Loeb I (f2007) (hereafter GL07) first noted 
that the tidal breakup of a stellar binary interacting with 
the SMBH can lead, at least for some orbital parame- 
ters, to a physical collision soon afterwards. If the two 
stars collide with a relative impact speed smaller then the 
escape speed from their surface, coalescence can occur, 
resulting in the formation of a new, heavier star^ 

Using Aarseth's direct integration scheme ([AarsethI 
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119991) to investigate the dynamics of the encounter, Gins- 
burg & Loeb found that the rate of colUsion events at 
the Galactic center is about ten times smaller than the 
rate of formation of hypervelocity stars (~ 10~^yr~^, 
iBrown et al.ll200 "6h. and that among the collision prod- 
ucts, stellar coalescence occurs in about twenty percent 
of the cases. In this paper we expand on their work and 
investigate the possibility that stellar encounters with the 
SMBH are responsible for the production of a population 
of rejuvenated stars around the SMBH. 

If a binary approaches the SMBH within its tidal dis- 
ruption radius, at the periapsis the stars can reach ve- 
locities a few percent of light speed, making relativistic 
effects potentially important in determining the prop- 
erties of the unbound population and the rate at which 
the stars are ejected at hypervelocities. Here, we perform 
the first A^-body simulations of binary-SMBH encounters 
that include post-Newtonian (PN) terms up to order 2.5, 
for both star-SMBH and star-star interactions. 

In §2 and §3, we describe the parameters we chose for 
our A'-body evolutions and their physical motivation. In 
§4, we discuss the likelihood of tidal capture of stars 
by the SMBH during close passages of binary systems. 
§5 and §6 consider in greater detail the populations of 
HVSs and stellar collision/merger products, respectively, 
including mechanisms like Kozai resonance that are crit- 
ical in understanding the long-term evolution of binary 
systems. A systematic comparison between Newtonian 
and PN integrations is provided throughout. Finally, in 
§6, we provide a summary and discuss future avenues of 
study. 

2. INITIAL MODELS AND NUMERICAL METHODS 

We use A^-body simulations to study the evolution of 
main sequence binaries as they make close passages by 
the SMBH at the Galactic center. We determine the con- 
ditions under which gravitational interactions produce 
HVSs, and the properties of the resulting distribution 
of both unbound and bound stars. All the simulations 
were carried out using the ARCHAIN (|Mikkola fc Merritt I 
120081) integrator, which includes PN corrections to all 
pairwise forces up to order PN2.5. The code employs an 
algorithmically regularized chain structure and the time- 
transformed leapfrog scheme to accurately trace the mo- 
tion of tight binaries with arbitrarily large mass ra tios 
(IMikkola fc Aarseth 112003 IMikkola fc Merritt ll200l . 

Given the significant computational resources per run 
made necessary by the high precision of our simulations 
we focus on a simplified set of initial conditions. We 
consider circular equal-mass binary systems with random 
initial orientations, initial separations ap — 0.05—0.2 AU, 
and individual stellar masses M* = 3 — 6 Mq. We define 
Mb to be the total mass of the binary system. Obviously, 
the choice of equal-mass components hmits our work. A 
treatment of binaries with a spe ctrum of masses simi- 
lar for instance to that found in iDuquennoy fc Mayor I 
()1991| ) is outside the scope of this paper, but could be 
included in a future work. 

We give the binary a tangential initial velocity Vi^ at 
a distance d — 0.01 — O.lpc from the SMBH of mass 
M, — Ax 10^ Mq, in effect setting the periapsis distance: 



for a Newtonian elliptical orbit. 

Our initial binary separations are close to the ex- 
tremes of the interval within which HVSs are produced. 
Since the ejection velocity increases with the internal 
binary energy, few stars will be ejected with velocities 
sufficiently large to escape from the Galaxy for oq > 
0.2AU, unless very larg e stellar masses are considered 
(|Gualandris et al.ll2005l ). We also note that, for the ini- 
tial distances d we choose, binaries with uq > 0.2AU are 
unlikely to survive for many orbits around the SMBH due 
to close encounters with field stars (Perots 2009a). The 
lower limit of 0.05AU is required by the consideration 
that smaller separations will result in contact binaries. 
GL07 pointed out that, assuming a constant probabil- 
ity per In(ao) for 0.02 < uq < 20AU, the probability 
to find a binary with a separation within the consid- 
ered interval is ~ 20%. We place the binaries on highly 
elliptical orbits about the S MBH, rather than e.g. on 
parabolic/hyperbohc orbits (| Gould fc Quillen 1120031 ). 

The origin of the young stars (and binaries) in the 
Galactic center is still a matter of ongoing research. One 
possibility is that stars/binaries formed in-situ, i.e. at 
distances of a few tens of milliparsecs from the SMBH. 
This model faces severe difficulties since tidal forces from 
the SMBH inhibit star formation at these distances. In- 
stead, it is typically assumed that the stars need to 
form further away from the SMBH and then migrate 
to the center within their lifetimes. Among the var- 
ious migration models that have been proposed, three 
produce orbits similar to those considered in this work 
(i.e. highl y eccentric): (i) the eccentric disk instability 
scenario f|Madigan et al. Il2008[ ) (ii) the cluster infall sce- 
nario aided by a n intermediate mass black hole (IMBH) 
(Gerhard "2001') and (iii) the triple disruption scenario 
(Perets 2009b). The first model considers a dynamical 
instability of an eccentric disk as the mechanism that 
drives the eccentricities of the stars away from their ini- 
tial values and produces near-radial orbits within a few 
Myr. The second model assumes stars and binaries form 
during the collapse of a giant molecular cloud and inspi- 
ral due to dynamical friction toward the Galactic cen- 
ter. The cluster, which is subject to tidal disruption, 
can reach distances of 10 — 50 mpc if it harbors an IMBH 
of mass ^ lO'^ Mq at its center. Tidally removed stars 
can then be scattered onto eccentric orbits by the IMBH 
([Merritt et al. Il2009l ). Finally, in the triple disruption 
scenario, young and runaway stars can form by disrup- 
tions of triples by the SMBH, which in some cases results 
in the capture of binaries on close orbits. 

At least some of the binary progenitors of observed 
HVSs were likely to have been part of the young stellar 
disk observed at the Galactic center, whoso inner edge 
lies a t a radius O.lpc from the SMBH (Lo ckmann et alTI 
|2008[ ). and this fact motivated our choice for the max- 
imum value of d. The distance 0.0 Ipc adopted as a 
minimum value for d is comparable to the radii of the 
smallest observed stellar orbits at the Gal actic center 
(|Eckart fc Genzel IIT997I : iGenzel et al. II1997I) . 

3. BASIC RELATIONS 
Three radii play a fundamental role in our simulations: 

• the periapsis separation Tpor of the initial orbit of 
the binary with respect to the SMBH; 
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• the tidal disruption radius rt of a single star due to 
the SMBH; 

• the tidal disruption (i.e. breakup) radius r^t of the 
binary due to the SMBH. 

A binary with initial separation ao will be broken apart 
by tidal forces if its center of mass approaches the 
massive object with in the distance ([Miller et al. .,20051 : 
ISesana et al. 1120091 ) 



rbt 



'Mb 



1/3 



ao- 



(2) 



As Hills (1988) noted, the tidal disruption of single 
stars is an important issue in the context of binaries. Be- 
cause of their low densities, main sequence (MS) stars do 
not survive close interactions with the SMBH; small ini- 
tial velocities take the stars too close to the SMBH and 
no HVSs are produced (compact objects, e.g., neutron 
stars and white dwarfs, can survive significantly closer 
approaches). Tidal disruption occurs for stars that ap- 
proach the SMBH more closely than rt, where 



(3) 



with i?, the stellar radius. We assume a mass-radius rela- 
tion R^/Rq = {M^/Mof-^^ ([Hansen et al. l[200l which 
yields i?* = 0.01 AU for = 3 Mq, and i?* = 0.016 AU 
for M^: = 6 M0. Stars with orbits meeting the criterion 
rnsr < n are fully disrupted (iLu minet & Carter 1986; 
lEvans fc K ochanek l ll98alFaber et al. II2005I) . Note that 
this condition is stronger than th e cla ssical Roche limi t 
formula given bv lPaczvnskTI ()1971f ) and lEggletoHl (|1983f ). 
Even when tidal disruption does not occur, we expect 
that mass will be stripped from the outer regions of any 
star passing within the Roche limit, which would have 
a fundamental impact on the subsequent evolution and 
orbit of the stars. 

For = 3Mq and M, = Mmw, n is ~ 1.10 AU, 
while for Af* = 6 M0 the corresponding value of rt is 
1.40 AU (i.e., MS stars with larger masses are disrupted 
at larger distances). 

A comparison of equations ([2]) and ([3]) reveals that for 
MS binaries we always have rbt > ft ( for contact binaries 
instead it is possible to disrupt a star before ejection 
occurs) . Our numerical simulations demonstrate that for 
large binary separations where rper < rbt, the two stars 
remain bound to each other up to distances rbt/2. 
This can be understood if the dynamical crossing time 
of the binary at rbt is of the order of the time scale of 
tidal breakup of the binary. Furthermore even at later 
times, after the binary is broken apart, both stars remain 
close to the original Keplerian orbit of the binary around 
the SMBH until reaching distances ~ rpcr. We will come 
back to this point below. 

Using equations ([T]) and ([3]) we can roughly estimate 
the lowest value of Win that avoids tidal disruption of the 
individual stars. The condition rpor > rt yields 



2GM, 



1/2 



(4) 



For d — 0.01(0.1)pc, this critical velocity is = 
44(4.0) km s-i for = 3M0 and 50(4.5) km s^^ for 
M* = 6M0. When PN corrections are included the or- 
bits are no longer Keplerian and the periastron distance 
can be found by including relativistic terms in the ex- 
pression for the potential (,Soffel II1989D : 



Vir) 



GM.Mh 



GM.L^ 



2Mbr2 c^M^r^ 



(5) 



where L is the binary angular momentum and c the speed 
of light. For a given initial velocity, the net relativistic 
effect consists of a decrease of rpor, yielding a higher limit 
for the initial velocity than for the Newtonian case. This 
value, in the case d = O.Olpc, is ^ 47 and 53kms~^ 
for M, = 3 and 6M0 respectively. For < 
would be tidally disrupted. 

We can derive an expression for the maximum initial 
velocity that should still be able to produce a HVS. The 
minimum condition is that the binary itself is broken 
apart. Together with equation ([2]), this gives 



vl stars 



Win < Wo 



2GM, 



1/2 



(6) 



(dVrbt + rf). 

Very few HVSs are likely to be ejected at > 

In this approximate model, HVSs are produced for 
Vi < Win < Between these limits, one member of 
a binary may be captured by the SMBH and begin to or- 
bit it at large eccentricity, while the other star is ejected 
with a velocity Wej that is larger than the escape velocity 
from the Ga laxy. An app roxii nate expression for Ve\ w as 
obtained bv lHillsl (fT988f ) and lYu fc Tremaine I (|2003l ): 



1770 



1/2 / Mu 



( ao \-i/2 / 
\0.1AU) \ 



2Mp 



1/3 



3.5 X lO^Mp, 



1/6 



/r km s" 



(7) 



with /r a function of the dimensionless closest approach 
parameter D = (rpor/ao)[2 M./IO*^ Mb]~^/^ An ex- 
pression for /r was derived bv lBromlev et al.l (|2006D : 

/r = 0.774 + (0.0204 + [-6.23 x 10"^ + {7.62 x 10"^ + 



(-4.24 X 10" 



.62 X IQ-^^ D)D}D]D)D, (8) 



which reproduces the spectrum of ejection velocities 
for binaries initially unbound relative to the SMBH. 
For bound orbits, the ejection speeds will typically be 
smaller, unless very large apoapses are adopted. 

4. TIDAL CAPTURES 

We first devised a set of simulations to determine the 
rate at which stars are captured by the SMBH. Here, 
"capture" means that a star passes either inside the stel- 
lar tidal disruption radius rt, or inside the innermost 
stable circular orbit (ISCO) at 3i?s, Rs = 2GM,/c^; 
RiSCO ~ 0.25 AU for M, = 4 x 10*^ Mq. Naively, one 
expects that passage of the binary within either radius 
would result in both stars being captured. In fact, due 
to the finite size of the binary, it is possible for one star 
to be destroyed and the other to escape as a HVS. 

We performed ~ 36000 integrations with equal mass 
M* = 3 M© binaries, oq = 0.05 or 0.2 AU and d = O.Olpc, 
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Fig. 1. — Probability that both stars (2P in the left panels) or only one star (IP in the right panels) of the binary fall within the tidal 
disruption radius rt (solid lines) and the SMBH ISCO SRg (dashed lines) as a function of the initial velocity The vertical dot-dashed 
line represents the value of iij. 



choosing random orientations and initial tangential ve- 
locities in the range 4kms~^ < Vin < 85kms~^. For 
each of the two binary separations, we considered both 
Newtonian and PN cases. Stars were treated as point 
masses, and we recorded all instances when a star fell 
inside 3i?s or rt. Given the possibility that only one of 
the two stars is captured during a close passage to the 
SMBH, it is important to include initial conditions with 
''per ^ ft (or equivalently Vm < Vi). 

Figure [T] shows the probability that one or both stars 
will pass inside 3i?s and/or rt as a function of the initial 
velocity Ui„. As a rule, passage of one star inside rt or 
3i?s implies that both stars are captured. 

The probability that both stars are tidally disrupted 
rises essentially to one at Vin ^ ~ 44kms~^ in the 
classical case, and at a slightly higher velocity in the 
relativistic case. 

The relativistic shift is well known: for a Keplerian or- 
bit with rper ~ 1 AU, the difference in the relativistic pe- 
riapsis with respect to the Newtonian value is ~ 0.12 AU, 
attributable to the attractive third term that appears in 
the relativistic potential in equation ([5|). Thus, when 
the binary has Win < 50kms^^, both stars are likely to 
b e tidally disrupted b y the SMBH. 

iGinsburg &: Loeb1 (|2006[ ) noted that the typical im- 
pact parameter that leads to the break-up of the binary 
by the black hole, is much larger than r^, the tidal disrup- 
tion radius for a single star. On this basis, they ignored 
stellar tidal disruptions. Our initial conditions are es- 
sentially the same as theirs; but as shown in Figure [T] 
these initial conditions would result in stellar tidal dis- 
ruption for a large fraction of the orbits, contrary to the 
assu mption of Ginsburg fc Lo eb. In their subsequent pa- 
per, (^niburg|&XoebJ((2003) assumed Vin < 25kms^^, 
which always leads to stellar disruption. Our study 
demonstrates that the relevant distance in the problem is 
the orbital periapsis of the center of mass of the binary 
(equation [Ij , which is also the distance of the closest 



approach of the stars to the SMBH. Indeed, for a wide 
range of initial velocities, the stars penetrate deeply the 
SMBH's potential weU. For Vin > 45kms"^ tidal dis- 
ruption is avoided but tidal perturbations of the SMBH 
would still be expected to have a significant influence on 
the subsequent evolution of the stars. 

Figure [1] also indicates that it is possible for only one 
star to be captured. In our calculations this mechanism 
does not produce a significant number of HVSs since this 
occurs for a narrow range of initial velocities: Win ~ Vi 
in the classical case, and Win ^ vi in the relativistic case. 
Between these cases we found that for wider binaries 
(oq = 0.2AU) the non-disrupted star is typically still 
on a bound orbit around the SMBH while for tighter bi- 
naries (oo = 0.05AU) the ejection probability is around 
^ 80% and the mean ejection velocity is ~ 3000 km s^^. 

If rpor is slightly larger than rt, the binary is usually 
disrupted at the first encounter with the SMBH. However 
we found a small fraction of orbits for which the binary 
survived for longer times. In our simulations, unbinding 
of the binary (i?b ~ 0) occurs approximately at 2 AU for 
ao « 0.05 AU and at ~ 10 AU for aa « 0.2 AU. In the 
subsequent evolution the separation between the stars 
does not change appreciably until the periapsis passage. 
After this point, the stars can either separate or become 
bound to each other again, continuing to orbit the SMBH 
as a binary. When this occurs, the minimum separation 
between the stars can become very small at later times. 

An example of the latter situation is illustrated in 
Figure [H which shows the evolution of the separation 
between the stars, as well as the orbit of the binary 
with respect to the SMBH, for an initial separation 
ao ^ 0.05 AU and initial velocity Vin = 51kms~^. After 
the first encounter with the SMBH, where the internal 
binary eccentricity suddenly changes, the periapsis sepa- 
ration between the stars is ~ 0.035 AU while the apoapsis 
separation is ~ 0.08 AU. After the second encounter, the 
stars orbit each other on an elliptical orbit with periap- 
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Fig. 2. — Example of an orbit whose periastron separation from the SMBH is Tper ^ rt- The binary is initially located at a: = 
2000AU(0.01pc) with a purely tangential initial velocity of = 51kms~^. The stellar masses are M« = 3M0 with initial binary 
separation ag = 0.05 AU. The left panel shows the orbit of the stars. The binary makes three revolutions around the SMBH before it is 
disrupted, after this point the individual stellar orbits are given by dashed lines. The upper inset panel schematically displays the different 
internal binary orbits during the first (A), second (B), and final loop (C). The lower inset panel shows the stars' orbits with respect to 
the SMBH during the first encounter. Note that the x and z-axes have very different scales. The right panel displays the change of the 
separation between the stars with time. The upper and lower limits of each "block" represent apoapsis and periapsis of the internal binary 
orbit respectively. 



sis separation ~ 0.028 AU. The third encounter with the 
massive object resuhs in a very eccentric orbit. FinaUy, 
after four encounters with the SMBH the binary is dis- 
rupted and the two stars begin to orbit the SMBH on 
two different high eccentricity orbits. 

The previous arguments suggest that the eventual 
shrinkage of the orbit after the encounter with the SMBH 
can easily trigger mass transfer and, in extreme cases, 
a coalescence between the stars. The merger remnant 
would orbit the SMBH on a highly eccentric orbit. In 
the case displayed in Figure [H the minimum separation 
between the stars is small enough to allow mass trans- 
fer after the second passage close to the SMBH. The 
Roche limit corresponding to the onset of mass transfer 
in binary systems of equal-mass stars would be ~ 2.8i?* 
while the dynamical instability limit in the case of co- 
rotating equal-mass MS s tars i s located around ~ 2.4i?* 

asio fc Shapiro I [1991 [T995h : for the case of eccen- 
tric binaries the radius of the Roche lobe is similar to 
that of circular bin aries when calculated at periastron 
(|Regos et al. 1 12005() . We also stress that this scenario 
would hold for higher values of rper, in which case it 
would not be restricted to small periastron separations, 
as will be shown below. 

From this first set of integrations we deduce that cases 
in which only a single star is captured are relatively rare, 
and occur when rper ~ ft- Stars ejected from these or- 
bits are expected to contribute weakly to the population 
of HVSs. This conclusion should hold in general when 
different values of d and M* are chosen. 

5. EJECTION OF HVS 

For the remainder of this paper, we focus on cases 
where Tpor > n. This allows us to neglect possible tidal 
effects on stars that would require a hydrodynamic treat- 
ment. 

The tidal disruption of binaries by SgrA* is generally 



TABLE 1 
Initial parameters. 



ao[AU] 


M.[M0] 


d[pc] 


Gravity 


0.2 


3 


0.01 


post-Newtonian and Newtonian 


0.2 


6 


0.01 


post-Newtonian and Newtonian 


0.05 


3 


0.01 


post-Newtonian and Newtonian 


0.05 


6 


0.01 


post-Newtonian and Newtonian 


0.2 


3 


0.1 


post-Newtonian 


0.2 


6 


0.1 


post-Newtonian 


0.05 


3 


0.1 


post-Newtonian 


0.05 


6 


0.1 


post-Newtonian 



accepted to be the main source of HVSs (stars ejected 
with V > 1000 km s^^) in our galaxy. The maximum 
velocity that can b e achieved with the classical binary 
supernova scenario aauw 1 IT961I) is < SOOkms-i for 
3 Mq stars. Interactions with a massive compact object 
seem necessary in order to explain the extreme velocities 
of HVSs. 

In this section we explore this mechanism by integrat- 
ing a set of ^ 10000 orbits featuring binary stars closely 
interacting with a SMBH. Two different stellar masses 
were considered: = SMq and 6Mq. The runs were 
performed for initial velocities that correspond to peri- 
apsis between rt and 6rbt • In all the simulations the final 
integration time was fixed to 60 orbital periods of the ini- 
tial binary orbit around the SMBH. For the adopted pe- 
riapsis distances, this integration time is long enough for 
the perturbations induced by the SMBH on the internal 
binary's orbit to grow dramatically. Table [T] summarizes 
the initial parameters we chose for the integrations. 

Because the binary is initially bound to the SMBH, 
no HVS is produced when two stars merge (i.e., if the 
relative velocity is lower than the escape velocity from the 
stellar surface after a collision) . On the other hand, if the 
stars collide without merging the ejection of a single star 
is still possible. Hence, in the following, we assume that 
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TABLE 2 

Ejection probabilities for the post-Newtonian (Newtonian) runs described in Table [T] for cases where rpcr < rbt- 



ao[AU] 


M*[M0] 


d[pc] 


Escape with 


Escape with 


Escape with 


Escape SMBH 








Dcj > 800kms-i 


> lOOOkms-i 


Dei > 1400kms-i 




0.2 


3 


0.01 


11.0(9.57) 


9.71 (7.35) 


6.60 (7.02) 


28.2 (24.3) 


0.2 


6 


0.01 


55.5(56.3) 


53.5 (49.7) 


27.7 (27.5) 


57.3 (60.7) 


0.05 


3 


0.01 


68.6(62.0) 


66.7 (62.0) 


65.6 (62.0) 


72.0 (62.0) 


0.05 


6 


0.01 


54.2(66.0) 


52.5 (65.1) 


52.5 (65.1) 


65.5 (66.3) 


0.2 


3 


0.1 


80.0 


74.5 


55.5 


89.1 


0.2 


6 


0.1 


80.9 


79.1 


71.8 


85.5 


0.05 


3 


0.1 


69.3 


69.3 


69.3 


71.6 


0.05 


6 


0.1 


44.3 


44.3 


44.3 


44.3 



a HVS ejection is possible even after the stars collide, 
unless the final product is a merger. 

5.1. Unbound population 

Figure [3] displays the average asymptotic ejection ve- 
locity of the stars unbound from the central object, as 
well as a comparison with the approximate model de- 
scribed by equations HJ, ([6]) and ([7|). Our results are 
in fair agreement with the predicted model, at least for 
d = O.lpc. When the apoapsis of the initial orbit is re- 
duced to d = O.Olpc and qq = 0.2AU, the ejection speed 
of the unbound stars can be a factor of 2 smaller than 
the typical values found for initially unbound stars. For 
qq = 0.05AU, on the other hand, the results are in good 
agreement with the model. The larger discrepancy in 
the case of wider binaries with respect to the theoreti- 
cal model (i.e., unbound orbits ) is due to their larger 
value of Tbt- Since the ejection velocity depends strongly 
on the Keplerian velocity of the binary when its compo- 
nents start to separate, the contribution of the binary 
orbital energy to Vej will be more important when the 
disruption occurs at larger distances from the SMBH. 
Another source of the observed discrepancy is due to the 
artificial truncation done at w^, beyond which a num- 
ber of unbound stars are still produced which represent 
the tail of the distribution with the lowest values of v^j ■ 
These cases are mainly due to binaries which are broken 
apart after few orbits around the central object. At each 
encounter the internal eccentricity of the binary changes, 
and there is a new chance for the stars to be separated 
and for a member to be ejected. 

To demonstrate the orbital evolution during an inte- 
gration. Figure 2] shows the change of the internal bi- 
nary eccentricity Cb as a function of the distance to 
the SMBH for several representative cases with differ- 
ent initial velocities in both FN and Newtonian grav- 
ity. The eccentricity was evaluated from the specific 
angular momentum and binding energy of the binary: 

e = -I- 2EL'^ /{GAhf. The initial separations, stel- 
lar masses, and apoapses were given respectively by 
ao = 0.05 AU, = 3 Mq and d = O.Olpc. The eccen- 
tricity behaves similarly for different initial conditions, 
increasing suddenly at ~ rbt, though in some cases the 
orbits may do several passages within rbt while maintain- 
ing a negative internal binary energy. When the binary 
has a small periapsis, several transitions between Cb > 1 
(unbound) and < 1 (bound) are possible. 

Results from our runs whose orbits have periapses 
within the binary tidal disruption radius (i.e., rpcr < ^bt) 
are described in Table [2] In general, we find for wider bi- 
naries with oo = 0.2AU, lighter binaries with = 3 Mq 



are less likely to produce HVSs capable of escaping the 
galaxy than heavier binaries with Ad^, = 6 M© , though 
they can produce a significant population of stars bound 
to the galaxy but unbound from the central SMBH. For 
tighter binaries with ao = 0.05AU the situation is re- 
versed, and lighter binaries are more likely to eject HVSs 
from the galaxy. For these tighter binaries, there are 
very few ejections of stars that unbind from the SMBH 
but remain bound to the galaxy, regardless of the stellar 
mass. 

A comparison between the velocity distributions of the 
ejected stars in the Newtonian and FN cases reveals that 
relativistic effects have a slight influence on the mean 
properties of the sample. In particular for ao = 0.05AU 
and d = O.Olpc the mean ejection speed is larger when 
the relativistic corrections are included. For larger sepa- 
rations, instead, this effect is absent and the simulations 
show similar results. However, we were not able to iden- 
tify a clear systematic effect due to the FN corrections 
in the properties of the ejected stars. 

5.2. The Bound Population 

Each star ejected during a binary-SMBH encounter is 
associated with a captured companion that loses energy 
in the process and becomes more tightly bound to the 
SMBH on a high eccentricity (e ^ 0.9) orbit. Bound stars 
are also produced after the tidal break up of a binary if 
neither star is ejected. The orbital parameters of the 
bound stars will be strongly correlated with the amount 
of energy carried away by the companion. For an initially 
unbound orbit the apoapsis distance of the bound star is 
approximately 

GM, 

a —, (9) 



and the corresponding orbital period is 



(10) 



The mean value of the semimajor axis A (with respect to 
the massive body) for three-body exchanges with equal- 
mass binaries can be expressed as (|Hills 1991) 



{A) w 0.56 



M. 

Mb 



2/3 



ao- 



(11) 



Bound stars can also be produced when the binary 
components merge. A coalescence remnant will not be 
able to escape the SMBH gravitational potential if the 
initial binary is bound to the SMBH unless significant 
mass loss occurs during coalescence. 
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Fig. 3. — Average ejection velocity of the stars unbound to the SMBH for stellar masses 3 Mq (filled circles) and 6 M0(open circles). 
The model described in §3 is superposed using dashed lines for 6 Mq and solid lines for 3 Mq. The horizontal dot-dashed line shows the 
approximate escape velocity from the Milky Way (~ 1000 km s~^). 



Figure [5] shows the dependence of the average final ec- 
centricity of bound stars on Vi^. A decrease in v-m (and 
thus an increase in the eccentricity of the initial orbit) 
results in a larger final eccentricity for the bound star, 
while larger stellar masses and initial binary separations 
decrease the eccentricity of the captured star's orbit. Fig- 
ure [S] displays the eccentricity instead as a function of 
the orbital period of the stars that remain bound to the 
SMBH at the end of our PN simulations (the Newto- 
nian calculations look similar). If we neglect the bound 
stars produced by mergers, three different families can 
be identified: 

1) Bound stars whose companion also orbits the 
SMBH (grey open circles in Figure [6]); 

2) Bound stars whose companion has been ejected by 
the SMBH but remains bound to the Galaxy (red 
triangles in the figure); 

3) Bound stars whose companion is a HVS escaping 
the Galaxy (filled circles in the figure). The ap- 
proximate escape velocity from the Galaxy was 
chosen to be 1000 km s~^ which corresponds to 
stars able to reach distan ces larger than ~ lOO kpc 
from the Galactic center (jKenvon et al. |[2008[ ). 

Due to the dependence of the apoapsis distance on the 
ejection velocity (see equation [9]) the three families oc- 
cupy distinct regions on the e — P plane. If more kinetic 



energy is deposited into the ejected star, the apoapsis 
distance of the bound star is reduced. In the case of 
equal mass binaries, the range of possible periods for the 
three populations of bound stars are independent of the 
stellar mass and the binary separation. 

On the other hand, in the case of unequal mass binaries 
with q — Afb/2Mi the apoapsis distance of the bound 
star can be approximated by 

GM, 

ai —q, (12) 

where Wej is given by equation [7| and for moderate mass 
ratios, the probability to be captured is almost indepen- 
dent by the stellar mass. 

Figure [6] shows that the stars with the highest values 
of P and e have reached velocities close to the escape 
velocity from the SMBH. Their companions also orbit 
the SMBH, but with much shorter periods. The region 
e — 0.92 — 0.95, strongly populated for ao — 0.2 AU and 
d — O.Olpc, is occupied by stars that were initially part 
of binaries with the largest values of Win. For d = O.lpc 
bound stars are mainly companions of ejected stars and 
the comparable region is empty. 

Stars whose the companion is escaping the Galaxy have 
orbital periods < 50yr. Many cases with orbital periods 
< lOyr are also observed. 

Finally, we note that in PN gravity the periapsis dis- 
tances of bound stars are about 0.1 AU smaller than in 
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log(R[AU]) 

Fig. 4. — Internal binary eccentricity as a function of the dis- 
tance R to the SMBH for different initial velocities (displayed in 
the panels). The binaries have initial separation ao = 0.05 AU, 
stellar masses SMq and the apoapsis is O.Olpc. Top: Newtonian 
integrations. Bottom: PN integrations. The vertical dot-dashed 
line gives the value of r^,^. For small the binary approaches 
very near the SMBH and the encounter results in the tidal break- 
up of the binary. In some cases the binary can complete several 
orbits around the SMBH before its stellar components are finally 
separated, even though rpci < fbt . Binaries with larger periapsis 
separations, or, equivalently, larger initial velocities, are likely to 
survive for a larger number of orbits. During each periapsis pas- 
sage, gravitational interactions produce significant changes in the 
internal eccentricity of the binary. 



the Newtonian regime. 

6. COLLISIONS AND MERGERS 

In this section, we study the evolution of binaries 
that are strongly perturbed during close passages to the 
SMBH, leading to stellar collisions and mergers. Stel- 
lar collisions due to either binary evolution or dynami- 
cal interactions are thought to be the main mechanism 
for the pr oduction of blue stragglers i n star clusters 
(ICollier et a l. 198i lBenz fc Hills Ill987t [Leonard lHOSOl : 
IMateo et al. Ill990f) . Similar processes may occur at the 
Galact ic center, producing a population of "rejuvenated" 
stars. iGenzel et al. I (|2003l ) proposed that the S-stars 



may be "super-blue stragglers" formed by collisions of 
lower mass stars and/or tidal heating. The appar ent nor- 
mality of their spectra fEisenah auer et al. 1120051 ) argues 
against rejuvenation ((Figer ,,20081 ) . but given the uncer- 
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Fig. 5. — eccentricity versus the initial velocity iiin of the initial 
orbit with respect to SMBH, for the orbits of bound stars with 
masses 3Mq (left panels) and 6Mq (right panels). The eccentric- 
ity is average over initial binary orientation. Higher stellar masses 
and lower binary separations tend to reduce the orbital eccentricity 
of the captured stars. 



tainties in the relaxation of merger products and the re- 
distribution of angular momentum, it is difhcult to draw 
firm conclusions as to their history. 

It is possible that the S-stars are drawn from the young 
population observed in the stellar disk that extends in- 
ward to within 0.1 pc of the SMBH, but while the S-stars 
are primarily B dwarfs, the young stars in the disk are 
mainl y luminous Wolf-Raye t and OB supergiants and gi- 
ants ([Paumard et al. Il2001f ). This limits the connection 
between the two stellar populations, as 0/W-R stars are 
typically more m assive and shorter-lived than B types. 
iLock mann et al.1 ([2008) proposed that binaries scattered 
from the stellar disk on highly eccentric bound orbits 
around the SMBH would be disrupted, ejecting one mem- 
ber as a HVS while leaving the other member bound. 
The inspiral of a cluster hosting an IMBH would gener- 
ate stars on comparable orbits. 

GL07 noted that the production of HVSs at the Galac- 
tic center could also result in collisions: the SMBH dis- 
rupts a binary, delivering an impulsive kick to one of the 
two stars at some direction with respect to the orbital 
plane, and under some circumstances the two stars will 
collide. We expect a dependence of the collision prob- 
ability on the initial binary orientation (unless gravita- 
tional focusing dominates) and a correlation between the 
number of HVSs produced and the number of collisions. 
However, GL07 concluded that the frequency of mergers 
is so low that no star in the Galactic center is expected 
to have been produced by binary coalescence. 

To investigate the rates of mergers and collisions us- 
ing iV-body techniques, we define the minimum impact 
parameter for a collision as 2R^ , and the two stars are as- 
sumed to merge when their relative speed upon collision 
is lower than the escape velocity from their surfaces. 

We find a much higher rate of coalescence than found 
by GL07, with most close collisions occurring only af- 
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Fig. 6. — Eccentricity versus period for the orbits of the stars that remain bound to the SMBH at the end of the PN integrations. The 
grey open circles correspond to stars whose companion is also orbiting the SMBH. Stars whose companion has been ejected by the SMBH 
but remains bound to the Galaxy are represented by red triangles, while stars whose companion escapes as a HVS arc indicated by filled 
circles. 



ter repeated encounters of the binary with the SMBH. 
The mechanism that leads to the majority of the close 
encounters in our integrations is Kozai oscillations. 

6.1. Kozai oscillations 

The perturbations on the inner binary's orbit caused 
by the gravitational interaction with the SMBH can re- 
sult in periodic oscillations (Kozai cycles) of both the in- 
te rnal binary e ccentricity Cb and the mutual inclination 
j ()Kozai Ill962f) . This occurs when the initial inclination 
Jin is sufficiently large and when the stars approach the 
SMBH at relatively large periapses ( ^ rbt) . More pre- 
cisely, the perturbations from the SMBH on the inner 
binary must always be weak and jin needs to satisfy the 
relation ic < jin < 180° — ic, where the critical angle i^ 
can be assumed to be ^ 40° in the case of initially cir- 
cular binaries and Newtonian gravity. Kozai cycles can 
result in a reduction of the periapsis separation, allow- 
ing the inner bodies to collide. The period of the cycles 
can be written in terms of the masses of the three bodies, 
the eccentricity of the outer binary e and their semimajor 
axis as: 

T = (1 - e ) I — I A(eb,Wb,Jin) 



Ph 



M, 



(13) 



where Ph is the periods of the inner binary , Pout the 
external period, Cb and cu^ the inner binary eccentricity 
and argument of periapsis, and K is generally of order 



unity l|Ford et al. Il2000n . Writing 
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the Kozai period becomes 
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Thus, in our integrations, the Kozai period is between ~ 
2 and '--^ 120 times longer than the orbital period around 
the SMBH, and the full effect of the Kozai cycles is only 
felt after several revolutions. 
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TABLE 3 

HVSs (ticj > 1000 kms"-"^), Collision and Merger frequency{%) for runs with rpcr < 6rbt- 



post-Newtonian (Newtonian); d = O.Olpc 




1 orbital period 


5 orbital periods 


60 orbital periods 




ao[AU], M4Mq] 


HVSs 


Coll 


Merg 


HVSs Coll Merg 


HVSs 


Coll 


Merg 


0.2, 3 
0.2, 6 
0.05, 3 
0.05, 6 


3.14 (3.00) 

15.2 (15.4) 

16.3 (16.8) 
15.6 (16.2) 


2.54(2.18) 
4.47(2.92) 
8.05(9.74) 
13.8(16.2) 


1.50(2.15) 
3.01(2.62) 
5.97(5.45) 
13.8(13.3) 


3.65 (3.25) 8.30 (9.10) 5.49 (5.29) 
16.5 (16.7) 18.3 (12.5) 9.11 (9.59) 
21.7 (18.2) 15.7 (23.4) 14.4 (20.6) 
15.9 (16.9) 24.3 (30.2) 23.9 (27.9) 


4.31 (3.26) 
17.4 (16.8) 
22.4 (21.0) 
20.8 (20.2) 


22.0 (20.5) 
30.6 (30.6) 
41.0(47.0) 
47.5 (55.5) 


17.0(16.6) 
27.5(27.4) 
40.1(43.8) 
47.1 (53.2) 


post-Newtonian; d = O.lpc 



1 orbital period 5 orbital periods 60 orbital periods 



ao[AU], M4Mq] 


HVSs 


Coll 


Merg 


HVSs 


Coll 


Merg 


HVSs 


Coll 


Merg 


0.2, 3 


29.0 


1.66 


1.00 


30.0 


5.00 


3.60 


30.0 


18.6 


17.1 


0.2, 6 


31.6 


4.33 


3.33 


35.6 


12.0 


9.33 


37.0 


20.3 


16.3 


0.05, 3 


25.0 


8.33 


8.05 


27.1 


21.7 


21.0 


28.0 


30.3 


29.6 


0.05, 6 


13.3 


18.3 


16.0 


14.3 


27.0 


24.3 


14.7 


40.0 


37.3 



If the orbit is initially circular, the maximum eccen- 
tricity achieved during a Kozai cycle is 

eb,max= [l-(5/3)cos2ji„]'/' , (14) 

the amplitude of the oscillations depending only on the 
initial inclination of the inner orbit relative to the ex- 
ternal perturber. The eccentricity of the outer binary 
remains roughly constant throughout the evolution, be- 
cause its variation is caused by octupole interactions, 
which are weaker than the quadrupole interactions driv- 
ing the oscillations of the inner binary's eccentricity 
dlakcda & Rasio 2005). 

Kozai resonances can be suppressed by additional 
sources of apsidal precession such as tides and relativis- 
tic precession in the inner binary. The inclusion of PN 
effects between the two stars damps the oscillations of 
the inner orbital eccentricity, eventually causing a de- 
crease in the n umber of obs erved collisio ns and mergers 
(|Holman et al. ,.1997 : Blaes et al. II2002D . 

6.2. Results 

The results of our runs show, for what we believe to be 
the first time, that binary break-up is not a prerequisite 
for physical collisions, and that the collision probability 
increases strongly with time: increasing the simulation 
time (and thus the number of revolutions around the 
SMBH) results in a dramatic increase in the number of 
collisions. The first point is of particular relevance for 
binaries with high values of Tpor, for which the SMBH 
more gradually reduces the internal binary angular mo- 
mentum. This would allow the two stars to closely in- 
teract through mass transfer and eventually merge while 
they are still bound to each other. We note that the first 
encounter with separation < 2i?* usually occurs just af- 
ter the closest approach of the binary to the SMBH, when 
the orbital parameters suddenly change (see Figure [J) . 
With regard to the second point: allowing the orbit to 
make only 2 — 3 extra revolutions around the SMBH in- 
creases the collision probability by up to three times with 
respect to shorter integrations (GL07). This is a conse- 
quence of the fact that the Kozai period is longer than 
the period of revolution around the SMBH. 

Figure [7] gives four representative examples of how the 
inner binary parameters change when repetitive collisions 



occur between the stars, with and without PN correc- 
tions. These cases are difficult to treat with an A^-body 
approach, but we may safely conclude that the stars are 
likely to merge after a short time. When relativistic 
terms are included, the maximum value of the binary 
eccentricity reached in a Kozai cycle, eb,max, decreases. 
However, we found that the effect of relativistic preces- 
sion on the likelihood of collisions and mergers observed 
in our simulations is weak. 

In Table [3] we list the fraction of observed collisions, 
mergers, and HVS ejections after the first binary-SMBH 
encounter, after 5 orbital periods, and after 60 orbital 
periods. After the first encounter, most of the binaries 
with small periapsis are broken apart, depositing one 
star on a tight orbit around the SMBH while ejecting 
the companion as a HVS. The frequency of HVS ejection 
is nearly unchanged over time because the orbits with 
larger periapsis do not get close enough to the SMBH to 
eject a member. Such orbits, instead, suffer strong per- 
turbations from the SMBH, altering the internal orbital 
parameters of the binary. The larger the periapsis, the 
larger the time required for the perturbations to become 
important, and the collision and merger frequency in- 
crease strongly with time. PN terms have important con- 
sequences only for initially small binary separations since 
the relativistic precession period is smallest for these sys- 
tems. However, as shown in Table |3] the fraction of col- 
hsions is reduced by no more than 10% for d = O.Olpc, 
since for tight binaries even small eccentricities lead to 
close stellar encounters. 

When the apoapsis is d = O.lpc the ratio between the 
timescale of relativistic precession in the binary and the 
period of the Kozai cycles becomes much smaller sup- 
pressing the Kozai resonance. Writing the timescale of 
relativistic precession as: 

"""af c2 (l-el) 

with c the speed of light, the product lot of the initial 
configuration gives the relative strength of relativistic 
precession to that of the tidal field of the SMBH. For 
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Fig. 7. — Evolution of the internal binary eccentricity and the mutual inclination j for different simulations in Newtonian (left panels) 
and PN (right panels) runs. In the cases shown, stars experience multiple collisions, and the time at which they finally merge is marked by 
a filled circle. For each periapsis distance considered , the corresponding Newtonian and PN integrations were started from the same initial 
conditions. The top three rows hav e d = O.Olpc and the bottom panels have d = 0.1 pc. The horizontal dashed lines marks the classical 
value fib njax given in equation II14I I. The vertical dot-dashed lines give the Kozai timescale. When the periapsis is sufficiently large the 
binary undergoes Kozai resonance, periodically changing both evj and j. When PN terms are included the amplitude of the oscillations 
decreases, resulting in a larger survival time for the binary and a smaller probability for a collision. The bottom panels show the evolution 
for systems with d = O.lpc, for which the stars do not merge if PN corrections are included. 



M = 6Mq and rpor > 12AU , this product can be as large 
as 10^ and the Koz ai cycles are completely suppressed 
(|Fabrvckv fc Treina ine 2007). The ejection probability, 
instead, increases due to the larger eccentricity of the ex- 
ternal orbit. The combination of these effects causes a 
smaller fraction of collisions/mergers respect to the num- 
ber of HVSs. 

Figure [H] displays the fraction of stellar mergers as a 
function of the periapsis distance rpor, with most cases 
occurring for Tpor > ^bt where the Kozai mechanism be- 
comes active. Figure [9] shows the cumulative fraction of 
stellar collisions, mergers and ejected HVSs in our PN 
integrations. Most of the HVSs are ejected during the 



first periapsis passage. Collisions and mergers, instead, 
typically occur later due to Kozai resonance, on a typical 
time scale t/2, where t is given by equation p3p . 

We note that strong mass transfer and/or merg- 
ers can lead to the formation of a rejuve nated star 
(|Vanbeveren et al. lITOOl IDrav fc Tout II2007D . 

The implications of these results for the actual rate 
of stellar mergers at the center of the Galaxy will de- 
pend on the dominant mechanism that drives binaries 
toward the SMBH. If binaries are scatt ered inward by 
"massive perturbers" (jPerets et al.ll200"7l ). most of their 
orbits will be weakly bound or unbound with respect to 
the SMBH, with apoapses of order parsecs, and will only 
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Fig. 8. — Fraction / of orbits leading to stellar mergers as a function of the periapsis distance for PN and Newtonian simulations. The 
vertical dot-dashed lines bound the region where HVSs axe produced. 
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Fig. 9. — Cumulative fraction of stellar collisions (solid lines), mergers (dashed lines) and ejected HVSs (dot-dashed lines) in our PN 
integrations. Most of the HVSs are ejected during the first encounter with the SMBH; collisions and mergers, on the other hand, occur 
more frequently later after the inner binary's parameters have evolved due to the gravitational perturbations induced by the SMBH. 
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encounter the SMBH once before being scattered onto 
different orbits. Given the low probability which we find 
for collisions during the first encounter with the SMBH 
for d = O.lpc, we do not expect a significant merger rate 
unless the binary population at the Galactic center is bi- 
ased toward low values of uq. Moreover, if the binary 
is initially unbound, the merger product is likely to be 
ejected at low velocities and it will neither remain bound 
(unless substantial mass loss occurs) nor become a HVS. 
Alternatively, mergers in this scenario may result from 
disruption of triples that can leave binaries on bound or- 
bits with small apoapse (|Perets 1120099 ). Obtaining an 
explicit value for the merger rate would be difficult in this 
case given the large uncertainties about the population 
of multiple stellar systems in the Galactic center. 

Alternatively, the binaries may form closer to the 
SMBH, perhaps in a disk, and migrate inward. In this 
case, the rate of collisions/mergers, as given in Table 3, 
would be comparable to the ejection rate of HVSs, i.e. 
~ 10~^yr~^. However, given the strong dependence of 
the collision probability on the distance of closest ap- 
proach to the SMBH, the merger rate will be strongly 
correlated with the periapse distribution of the infalling 
binaries. The generation of massive new stars through 
binary-SMBH encounters in this scenario could possi- 
bly explain the presence of the stars that appear to be 
the youngest a mong the S-sta r gr oup at the Galactic 
center (see also iPeretsI (|2009b[ ) and iPerets fc Fabrvckv I 
([2009) for discussion on rejuvanting stars through Kozai- 
oscillation induced mergers). 

In our simulations, the distribution of the semima- 
jor axes and eccentricities of the captured stars are 
fixed by the initial distance d of the binary we adopt, 
and therefore is completely arbitrary. Nevertheless, as 
previously discussed, the apoapses found in our runs 
are consistent with those predicted assuming the binary 
progenitors of HVSs were scattered on eccentric orbits 
bound to the SMBH. The orbital eccentricities we find 
(e > 0.95), it should be noted, are larger than those ob- 
served for the S-stars at the Galactic center, which are 
consistent with a th ermal distribution (i.e., N{e) oc e). 
IPerets et al. I ()2009l ) have shown how perturbations from 
a dense cluster of 10 M© stellar black holes, expected to 
accumulate around the SMB H due to mass segregation 
(jHopman fc Alexander |[2006[ ). tend to randomize orbits 
with high eccentricity in ~ 20Myr. Alternatively, if 
the binaries are transported to the center by the inspiral 
of a massive cluster containing an IMBH, interactions 
with the sam e IMBH will therma lize the eccentricities 
in IMyr (|Merritt et al. II2009D . These processes of 
orbital randomization result also in the evolution of the 
periapsis distance rpor, but on a timescale much larger 
than that required for the binary systems to collide. We 
note, finally, that the presence of external forces on the 
inner binary as well as tidal friction between the stars 
could suppress the Kozai resonances reducing the num- 
ber of stellar collisions. 

6.3. Effect of the Initial Binary Orientation 

Based on the results of our simulations, we find no 
strong dependence of the collision and/or merger proba- 
bility during the first SMBH passage on the initial binary 
orientation. To quantify our results, we use the Rayleigh 
(dipole) statistic Cc (Cm) defined as the length of the re- 



sultant of the unit vectors h, i = 1,...,7V, where h is 
perpendicular to the initial internal orbital plane of the 
jth |3jjja,ry system whose component stars collide (merge) 
during the integration and N is the number of collisions 
or mergers. Here the subscripts c and m refer to 'colli- 
sion' and 'merger' respectively. 

For an isotropic distribution C ^ VTv", while C ~ N if 
some preferential direction exists (|Ravleigh Ill919[ ). 

In order to test the agreement with an isotropic distri- 
bution we used a Monte Carlo approach: for each set of 
simulation parameters we generated a set of 1000 sam- 
ples, each containing N randomly oriented orbits. Then, 
for the whole sample we computed the average of C and 
the corresponding standard deviation and we compared 
the resulting values with those obtained for the collid- 
ing/merging binaries. 

Table |1] gives the estimated values of C after the first 
binary-SMBH encounter as well as the expected values 
for an isotropic distribution. Assuming a confidence in- 
terval of 90%, most of the values shown in Table [4] are 
consistent with an isotropic distribution. 

We find values of Cm too large to be consistent with 
isotropy only for PN integrations where qq — 0.05 AU 
and d = O.Olpc, though for Newtonian integrations with 
the same apoapsis and oq = 0.05 AU , the values of Cm 
are quite large as well, possibly indicating some degree 
of anisotropy. In Figure [101 we plot the spatial projec- 
tion of the angular momentum eigenvectors for the in- 
ner binary, showing that when the angular momentum is 
aligned with the -l-y-axis, the stars are unlikely to merge. 
These binaries are retrograde and their internal orbital 
plane coincides the orbital plane of the binary with re- 
spect to the SMBH. In all the other cases we conclude 
that the collision probability does not show any signifi- 
cant dependence on the initial binary orientation during 
the first encounter with the SMBH. 

At later times stellar collisions are mainly produced 
by the Kozai mechanism which is active only for large 
initial inclinations jin . This corresponds to a low prob- 
ability of collisions for binaries with the angular momen- 
tum aligned with either the iy-axis. For these cases 
our Rayleigh statistic results are still consistent with 
isotropy, but this is found to be an artifact attributable to 
the dipole nature of the Rayleigh statistic. A more care- 
ful analysis of the orbital distribution of binaries whose 
components merge at later times shows that, as expected, 
the population of colliding binaries is strongly biased to- 
ward large inclinations of the inner orbit with respect 
to the external orbital plane. An illustrative example is 
given in fig llll which shows the projection of the internal 
angular momentum vectors for binaries with gq — 0.2 AU 
and d — O.Olpc merging between 10 and 60 orbital pe- 
riods. Even though the result of our statistical analysis 
in this case is consistent with isotropy, the distribution 
of the angular momenta is clearly anisotropic with most 
of the collisions occurring for large inclinations. The dis- 
crepancy is a consequence of the dipole nature of the 
Rayleigh statistic: the majority of the unit vectors U 
have a quasi-antialigned counterpart giving a net contri- 
bution close to zero. 

7. SUMMARY AND DISCUSSION 

We have carried out post-Newtonian numerical inte- 
grations of binary stars on highly elliptical orbits around 
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TABLE 4 

Raleigh statistics Cc and (^m and values expected for an isotropic dihtribution for collisions (EVc) and mergers(EV^) 

AFTER the first DINARY-SMBH ENCOUNTER. 





Post-Ncwtonian(Newtonian) ; d=0.01pc 




an f AC/1, M,[M«1 


t IN EV 
sc./ ^ y c 


t IN 


EVm 


0.2, 3 


0.33 (0.30) 0.18±0.15 


0.25 (0.29) 


0.23 ± 0.18 


0.2, 6 


0.19 (0.23) 0.17±0.14 


0.21 (0.25) 


0.20 ± 0.16 


0.05, 3 


0.20 (0.18) 0.12±0.10 


0.26 (0.21) 
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0.05, 6 
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Fig. 10. — Projection of tiic angular momentum eigenvectors for binaries with ao = 0.05 AU and d = O.Olpc merging during tiie first 
encounter witii the SMBH for PN (left panels) and Newtonian (right panels) runs. A slight dependence of the merger likelihood on the 
initial binary orientation is observed, with a deficit apparent for initial orbits whose angular momentum is aligned with the +i/-axis. The 
deficit is less significant in the Newtonian integrations. 



a supermassive black hole and classified the final fates of 
both stars. Our main conclusions can be summarized as 

follows: 

1. If the binary orbit passes through the stellar tidal 
disruption radius rt of the SMBH, both stars are likely 
to be disrupted during the initial passage, although there 
is a non-trivial chance that one star may remain outside 
the tidal disruption radius while the other passes within. 
We found that after the binary internal energy becomes 
positive the stellar orbits remain close to the initial tra- 
jectory of the disrupted progenitor. This demonstrates 
that HVSs can suffer strong tidal perturbations and mass 
loss before ejection. 

2. Binaries with pcriapscs in the range rt < rpcr < rbt 
with ~ rbt the tidal breakup radius of the binary are 
the main hypervelocity star progenitors. In these cases 



the (hypervelocity) ejection probability is typically larger 
than 50%, although the combination of wider separations 
(oo = 0.2AU) and small apoapscs id = O.Olpc) can result 
in probabilities less than 30%. 

3. For tighter binaries (taken here to have separations 
of 0.05 AU), there is a strong inverse correlation between 
the pcriapsc distance of the binary and the HVS ejection 
velocity, while the dependence is present but weaker for 
wider binaries. 

4. While many binaries perform several orbits around 
the SMBH before merging and/or colliding, the vast ma- 
jority of HVSs were produced during the first passage 
by the SMBH. This is due to the strong dependence of 
HVS production on the pcriapsc distance of the binary 
orbit: binaries that get sufficiently close to the SMBH 
are broken apart, while those farther away will typically 
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Fig. 11. — Projection of the angular momentum eigenvectors for 
binaries with oq = 0.2 AU and d = O.Olpc merging between 10 
and 60 orbital periods for Newtonian runs. This shows how the 
Kozai mechanism selects highly inclined orbits (the external orbital 
plane lies on the x-z plane and the direction of angular momentum 
corresponding to the external orbit is given by a blue continue 
arrow) . 

be unbound or collide. 

5. The tidal break-up of equal-mass binaries by the 
SMBH results in a population of bound stars with peri- 
ods and eccentricities defined mainly by the fate of the 
companion star (Fig. 6). Companions to escaping hy- 
pervelocity stars remain bound to the SMBH with or- 
bital periods P < 50yr, but also extremely short periods 
P < lyr are found. A much larger range of periods 
(10^ < P < 10^) is instead possible for stars whose com- 
panion is also left orbiting the SBH. 

6. For sufficiently inclined orbits, the Kozai mecha- 
nism produces periodic oscillations of the inner binary 
eccentricity. This reduces the internal binary periap- 
sis allowing for close stellar encounters. We find that 
the probability for collisions between the components of 



the binary increases with time, resulting in substantially 
larger numbers of mergers than in simulations that stop 
after the first passage. 

7. When the external orbital period is large compared 
to the period of relativistic precession in the inner bi- 
nary the PN corrections tend to suppress the oscillations 
reducing the probability of a collision which, however, re- 
mains high in our integrations. For initially bound orbits 
with Tper < 6rbt the probability that the two stars collide 
is similar (for d = O.lpc ) or much larger (for d — O.Olpc) 
than that of hypervelocity ejection. 



Stars caught on highly eccentric orbits around the 
SMBH (with periapsis rt < r < 2rt) can eventually 
spiral in as their orbital energy is converted into heat 
thro ugh tidal interactions during each periapse passage 
(e.g. [Alexander &: Morris l[20bS) . Such stars would have 
orbital properties generally consistent with the S-star 
population, but A'^-body calculations cannot accurately 
predict the properties of the merger products. Instead, 
in a future work, we plan to study the evolution of bi- 
naries as they orbit the BH and collide or merge using 
a smoothed particle hydrodynamics (SPH) code, to de- 
termine the properties of the resulting stars, including 
their masses and angular momentum distributions. This 
should help to clarify whether they would be expected 
to possess normal spectra or features more commonly 
associated with blue stragglers. 

Finally, our study provides concrete evidence that rel- 
ativistic effects can have important consequences for the 
distribution of HVSs and the collision/merger probabili- 
ties of HVS progenitors. 
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